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Abstract 


An examination of the haze regime shows that: (1) the aerosol properties of 
a surface meteorological range generally affect a mixing layer to 5 km altitude, 
and (2) the lower and upper visibility limits of a haze regime are defined by me- 
teorological ranges 1,2 km and 15 km respectively, Within these limits eight me- 
teorological ranges are selected for developing uv, visible, and ir aerosol attenua- 
tion coefficients. An aerosol scale height is derived for each meteorological 
range, Finally, the computed aerosol attenuation coefficients are presented as 
tabulations, which include previously published attenuation parameters (aerosols, 
molecules and ozone) to 50 km altitude, 
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Symbols 


Horizontal path length (km) 
Aerosol scale height (km) 
Altitude (km) 

Aerosol index of refraction 
Aerosol number density (em73) 

Constant proportional to total number of particles between ry mie To 
Molecular number density (cm) 

Meteorological range (km) 

Particle radius (microns) 

Horizontal transmission 

Transmission between sea level and altitude h 

Transmission between altitude h and space 

Transmission between two altitudes above sea level 

Atmospheric ozone absorption coefficient (km7}) 


Aerosol attenuation coefficient (km7}; 
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Rayleigh (molecular) attenuation coefficient (km) 


Extinction coefficient (km?) 

Zenith angle 

Wavelength (microns) 

Aerosol scattering cross section (cm?) 

Rayleigh scattering cross section (em?) 

Ozone optical thickness from sea level to altitude h (0-h) 
Ozone opticai thickness from altitude h to space (h-o) 
Aerosol optical thickness from sea level to altitude h (0~h) 
Aerosol optical thickness from altitude h to space (h-~) 
Rayleigh optical thickness from sea level to altitude h (0-h) 


Rayleigh optical thickness from altitude h to space (h--») 
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Extiretion optical thickness (molecular + ozone + aerosol) from sea level 
to altitude 4 (0-h) 


Extinction optical thickness (molecular + ozone + aerosol) from altitude h 
to infinity (h-°) 


y Aerosol size distribution friction 


Vertical-Attenuation Model With Eight Surface 
Meteorological Ranges 2 to 13 Kilometers 


1, INTRODUCTION 


A series of atmospheric attenuation parameters which vary with wavelength 
and altitude are useful for carrying out a variety of exploratory calculations, Such 
information can take the form of curves, tabulations, or analytic expressions. It is 
recognized, however, that limitations exist due to variability of the atmosphere's 
constituents, expecially the aerosol content of the lower troposhere, which contrib- 
utes extensively to the optical thickness. For example, in the photopic region, as- 
suming a representative wavelength A =0.55 and a meteorological range of 23 km 
near the surface, the aerosol content in the first 3 km above sea level accounts for 
about 70 percent of the total optical thickness. If surface conditions are hazy or 
polluted, the aerosol content accounts for a larger nercentage. This suggests that 
the treatment of atmospheric attenuation can be improved by introducing aerosol 
parameters related to the easily measured meteorological range, that is, by in- 
troducing quantitatively a haze regime which is considered as encompassing a 
series of meteorological ranges between those associated with normally clear con- 
ditions and fog, See Figure 1. ; 

Near the surface and at low altitudes, the aerosol constituent is ubiquitous and 
highly variable, and an aerosol component is present even for a very clear atmo- 
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aspheric condition, Thus, the boundaries of the haze vegime (used in the sende of 
diminished meteorological range) muat be defined, The hage limita will be bawed 
on the Koschinieder (1924) definition, 
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ext 


where 
7 (1b) 


Vy is the meteorological range (km) and Boxt’ Bie By are the extinction, Rayleigh, 
and aerosol attenuation coefficients (km~!) respectively, Equations (1a) and (1b) 
apply to sea level conditions and the photopic region represented by A = 0.65 y, 
When converted to a log-log trace (Figure 1), several boundaries can be designated 
conveniently. The Rayleigh limit, 336 km, is based on the standard atmosphere 
(15°C, 1013 mb), Foitzik (1938) found that the haze-fog transition is relatively 
abrupt, and, therefore, can be readily identified; and that Vy * 1,2 km represents the 
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(rataition, Pollak'a peaall ta vonfticmed adequately hy Netburger and Chien (1nd0), 
Mualtedeh (ENN), Clow (UNA, PUOd and the cnalyada by Midmidge (1adn), 

Trvonteaal, (he Uterature pestaining to a boundary condition between a clear 
ant hagy atmonphere ia meager, Since no abrupt changes voor in thia region, the 
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The meteorologival vangea and corveaponding parametera shown in Migure | and 
Table 1 will be waed in the material to follow because they are apaced at convenient 


Table 1, Meteorelogical Ranges and Corresponding Parametera 
(Representative Photapic Wavelength A = 0.88) 


Vy Poxt, a a} Ph ad Hp 

(Km) {km *) (km ") (km) (km) 
P} 1,988 O,0116 1.943 0,84 
N| 1,303 0.0116 1,201 0,80 
4 0,078 0.0116 0,966 0,98 
§ 0,782 0,0116 0,770 0,99 
8 0,682 0.0116 0,640 1,03 
a 0,489 0.0116 O.4TA 1,10 
10 0,391 0,0116 0,379 1,15 
13 0,301 0,0116 0,280 1,23 

Vy - meteorological range 

Pon - extinction coefficient 

B. - Rayleigh attenuation coefficient 

By - aerosol attenuation coefficient 

H - aerosol scale height 
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logarithmic intervals, are adquately separated from the haze-fog transition, and 
are within the hase regime characterised by diminished meteorological range, The 
scale height in the last column of the tabulation will be discussed later, 


2 SPROTRAL METEOROLOGICAL RANGES 


The concept of photopic meteorological range can be widened spectrally if con- 
current measurements of attenuation coefficients at other wavelengths are avail- 
able. The work of Curcio, Kneastrick, and Cosden (1961), which is the basis for 
Figure 2, is an example where 8 (V4. A) was obtained through a series of concurrent 
measurementa, Because of the quantity of data obtained, the results for the me- 
teorological range \y = 4 km are considered representative by these authors. A 
family of distributions, 6, (V, »A), can be computed if the B, (V4 , A) values are used 
in conjunction with Eq. ()) so that 


3.91 3.91 
BV) Ad = BA(Vgs A) ks - 8, a0] / [A F059 | ' (2) 
n 4 


V,, (km) being the photopic (A = 0.55) meteorological ranges of interest. Using 
kq. (2), the aerosol attenuation coeffic.ent is found for various combinations of me- 
teorological range and wavelength, 0.27 to 2.17y (Table 2), that is, 160 surface 
values, 6 (VY, ,A). The shapes of the distributions so determined (Figure 2) con- 
form rigorously to the distribution for V_ = 4 km, on which they are based. 

Because of the functional importance of Eq, (2), it would be in order to examine 
its implications, especially those related to particle size considerations. If we 
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Figure 2. Wavelength Distributions of 
the Aerosol Attenuation Coefficient for 
= 6 km and Vy = 10 km Derived From 
Vy = 4 km Using Eq. (2). The Vy = 4km 
curve is obtained from measurements 
by Curcio, Knestrick and Cosden (1961), 
which included the wavelength region 
0.40s A< 2.174. An extrapolation to 
0.27u permits computations for an 
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consider a real atmosphere the aerosol sizes within unit volume determine the 
aerosol attenuation coefficient described by 


Table 2. Surface Aerosol Attenuation Coefficients Corresponding to Figure 2 
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" 
Bp (m, A) = f ¢, (m,r,A) n (r) dr, (3) 
ry 
n(r) * Ny (V,,) Wy (r) (4) 


and when combined 


by) 
2 

BB (r,A,V,,) = ON. (V,)) [ % (v,A) W(r) dr. (5) 
1 i 


in these expressions, 6. is the aerosol attenuation coefficient; the index of refrac- 
tion is m (to be omitted following Eq. (3) because subsequent considerations will 
assume m invariable); ry and Tp are the lower and upper radii limits of the size 
distribution n(r); N a is a constant proportional to the total number of particles be- 
tween ry and Poi yw is the size distribution function (the same for all selected me- 
teorological ranges), 

It is implicit in Eq. (5) that B p and corresvondingly the aerosol number density 
determine the meteorological range (V,)- The integral in Eq. (5) is a wavelength 
function indevendent of the meteorological range, which accounts for the conformity 
in shape of the curves in Figure 2, 


3. STATEMENT OF OBJECTIVES 


The material thus far has dealt with: (1) the limits of the haze regime (in 
terms of meteorological range), (2) derivation of spectral aerosol attenuation co- 
efficients (0.27 ¢ A< 2.17) for a series of meteorological ranges, and (3) an ex- 
amination of the assumptions implicit in the derivation of these coefficients. Now, 
a statement of objectives can be made. Specifically, aerosol scale heights will be 
determined for the coefficients in accordance with their meteorological range and 
their vertical distribution. Then values of the coefficients for km intervals (0-5 km) 
will be computed. To the coefficients will be added previously published Rayleigh, 
ozone, and aerosol parameters for altitudes to 50 km (Elterman, 1968) in order to 
formulate an attenuation model for a haze regime with eight meteorological ranges 
(2< Vs 13 km). 


4. AEROSOL MIXING LAYER 


The procedure for assessing the aerosol attenuation coefficient above the surface 
can parallel that used for a clear atmosphere (Elterman, 1968), which entaila the ap- 


t 


ip 
AS 


7 


plication ofa suitable aerosol scale height. As an introductionto scale height con- 
siderations, it is noted that, meteorologically, a significant role is assigned to the 
altitude interval up to several km above the surface, This is a region of strong 
vertical mixing determined by such factors as heat-transfer across the earth-air 
j interface, advective winds, and consequent turbulence attributable to the region's 
topography. The resultant vertical flow, mechanical and convective, is character- 
ized, meteorologically, as a mixing depth equivalent to the vertical extent of the 
mixing layer. When dealing with aerosol attenuation coefficients, aerosol condi- 
tions in this layer can be examined in terms of mixing ratios such that for a se- 
lected altitude h 


By (h, A) a, (A) Ny (h) ih 


6. (h, A) o., (A) N. (h) 
where B, and B.. are respectively the aerosol and Rayleigh attenuation coefficients 
(om7}); and Ny and N,. are respectively the aerosol and molecular number densi- 
ties (em), The terms @. and O.. which are respectively the aerosol and Rayleigh 
cross sections (cm), tend to remain constant with altitude (a reasonable assump- 
tion), Eq. (6) then asserts that B.(h,A) / B. (h,A), known as the optical rnixing ratio, 
is proportional to the number density mixing ratio Ny (hh) /N, (h). The size distri- 
bution for N_ comprises aerosols sufficiently small to be responsive to the usual 
factors conducive to mixing. Meteorologically, the mixing depth is considered to 
be 3 km or less, so that an aerosol mixing depth determined from an optical mixing 
ratio or a number density mixing ratio or even by direct (in-situ) measurement of 
N_(h) should be in agreement, 
The conclusion, however, based on aerosol measurements sufficient to provide 
a meaningful average, is that the mixing depth normally extends to a greater alti- 
tude, Siedentopf's (1944) sky luminance measurements (18 aircraft flights) show 
that on the average, the aerosol concentration decreases exponentially with altitude 
| and that the scale height undergoes a significant change between 5 and 6 km. 
Penndorf's (1954) analysis of soiar attenuation observations (8 aircraft flights) 
| shows the scale height change to occur at 4.5 km, An examination of Rosen's (1967) 
: balloon photoelectric countermeasurements, selecting only those profiles where 
the surface layer is readily discerned (37 profiles obtained on ascent, descent and 
for 2 wavelengths), shows that the average mixing depth is 5.4 km, An analysis of 
optical probing measurements (Elterman, Wexler, and Chang, 1969) yielded optical 
mixing ratios which show that the depth of the surface layer averages 5.3 km (79 
profiles at 0.55u wavelength). Blifford and Renger (1969) completed a series of 22 
eerosol collections using an aircraft-mounted impactor. Samples obtained to 9.1 km 
provide mixing ratios that indicate the mixing depth to be in the rangr 3-6 km, 


Measurement of the atmospheric aerosol distribution has received considerable 
emphasis in the USSR, for example, in the work of Faraponova (1965), who conducted 
more than 60 aircraft flights in cloudless weather (solar atmospheric attenuation) 
to 6.5 km altitude, and in the summary by Kondratiev (1969). In general, the USSR 
findings are compatible with those previously described. 

An overall assessment of the results shows that the aerosol content for the Jow 
altitudes is characterized by a mixing depth between 4,5 and 5.5 km, As has been 
mentioned, it is somewhat higher than the mixing depth of the meteorology disci- 
pline. However, the difference is understandable when it is considered that in 
almost all instances, aerosol measurements were conducted over the continent and 
with cloudless skies, whereas the mixing layer in the meteorological sense repre- 
sents less limited conditions, Accordingly, the designation ''aerosol mixing layer" 
will be used, and will be assigned a depth of 5 km (considered representative), 
Within the aerosol mixing layer, considerable variation and stratification (frequently 
due to inversions) can occur but, on the average, the particle distribution, N,(h), 
and correspondingly the aerosol attenuation coefficient, B (h), decreases exponen- 
tially for the altitude region 0-5 km, The rate of decrease can be expressed as a 
constant scale height although, as will be shown, not necessarily the same scale 
height for each meteorological range. At 5 km, effects of mixing are substantially 
diminished. 


5. THE USE AND JUSTIFICATION FOR A SCALE HEIGHT FAMILY 


If aerosol conditions at the upper terminus of the aerosol mixing layer are rel- 
atively stable compared to those at lower altitudes, as discussed, it suggests that 
the scale height characterizing the aerosol mixing layer is related to the meteoro- 
logical range. The existence of such a relationship was examined by means of 
Figure 3, using A = 0,55. Specifically, the scale heights were determined by: 

(1) utilizing the surface values, By (A 55° Vy) for the meteorological ranges of in- 
terest derived from Eq. (2); and (2) taking from published tabulations (Elterman, 
1968) at 5 km the aerosol attenuation coefficient, B, (fs A 55) = 5,0 X 1072km7?, 
This quantity (assumed representative because the tabulations are based on 79 sets 
of measurements) is considered relatively independent of the meteorological range 
for reasons already given, With surface valves and the 5 km value established, the 
aerosol scale height (H,) was derived for each meteorological range by using 
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Figure 3. Relationships of Four Aerosol Scale Heights q 
(Hp) With Meteorological Ranges (Vy), Aerosol Attenua- 

ticn Coefficients (Bp) and Aerosol Mixing Layer Altitudes 

(0-5 km). The aerosol scale height family was computed 

using A = 0.554. The dash line represents values of 6 pi) 

above 5 km (Elterman, 1968) 


It is implicit in this procedure that the derived scale height family is independent 
of wavelength, With Ho known, the values for B (bd) can be calculated (using 

Eq. (7)) for each km interval up to 5 km (the aerosol mixing layer), To summarize, 
the decrease of B_ with altitude is represented by a family of scale heights and each 
scale height depends on the meteorological range of interest (Table 1), 

The validity of this procedure depends not only on the value but also on the 
variability of Bi (hs. r) as it affects the related parameters of interest. In this 
respect, we note that the aerosol optical thickness up to 5 km, 1 (hos, A), is an 
important objective; and also that it is obtained by integration, which makes it 


sensitive to change of its composite elements, especially at low altitudes. Thus, a 
suitable evaluation of the aerosol scale height family derived in accordance with 
Figure 3, is to vary By (hs, av 55) significantly and to examine the optical thickness 
and transmission change for the first 5 km and for the several meteorological 
ranges indicated. This was done by changing 6 p(hg+A : 55) = 5.0 X 10°3km 7} by a 
standard deviation, 0= + 3.4 107%, For this melculation, as previously mentioned, 
the mean By and o values at 5 km altitude were obtained from 79 selected optical 
probing measurements (Elterman, 1968), The resulting aerosol optical thickness 
and transmission changes are shown in Figure 4, For the eight meteorological 
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Figure 4. Comparison of Vertical Opti- 
cal Thickness and Transmission for +0 
(Standard Deviation) and a Mean Aerosol 
Attenuation Coefficient Bp (h5, 4,55) at 
the Top of an Aerosol Mixing Layer Hav- 
ing 5 km Depth, T and 7 represent the 
transmission and optical thickness based 
on Bp (hg, A,55) = 5.0 X 10-Skm"+, the 
mean of 79 measurements. T, and Ty 
represent the transmission and optical 
thickness based on 6, (hs, A,55) - 9% 
where 0= 3,4 X 10-3," To and To repre- 
sent the trans:nission and optical thick- 
ness based on Bp(hs, A,55) + 0 
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ranges, a +0 of the aerosol attenuation coefficient corresponds to an average of 

3.6 percent change in vertical transmission (attributable to aerosols only), Simi- 
larly, a ~O corresponds to a 6.4 percent verticai transmission change. The changes 
are unequal due to exponential relationships. These changes are relatively modest 
especially when it is considered that the aerosol attenuation coefficients represent- 
ing the meteorological ranges 2 km to 13 km extend over an order of magnitude, 
Based on this discussion, it is concluded that for conditions of diminished meteoro- 
logical range in the aerosol mixing layer, there is sufficient justification for the 
use of an aerosol scale height family, and that the scale height selected is best de- 
termined from the meteorological range of interest, 


6. THE OPTICAL THICKNESS COMPUTATION 


Since the aerosol attenuation coefficients, B, (ALA, V, ), as well as the aerosol 
scale heights, Hy (constant for each meteorological range) are known, an analytical 
expression for the aerosol uptical thickness, Tp! was derived as follows; 


h 


. ; 8 
Tp (HA, Vp) f B,(n,A,V,) dh (8) 


For a given wavelength, the aerosol scale height expression applicable to the haze 
regime depends on the meteorological range and the vertical distribution of the 
aerosol attenuation coefficient, Accordingly, 


peep ae Srmemportltg 7 oe af ea i d , 


-h/ HOY,) 
= . h - h . » (10 
Tp (hd, V,) Hy (V,) B, 9° Vy) H, (V,) B,( 5 A,V,) e (10) 
Applying Eq, (9) to the bracketed factor in Eq. (10), 


T(t AV,) = Hy (V,) (8, (h.4V,) - BAY.) (11) 


Equation (11) was used to compute the aerosol optical thickness for the combina- 
tions of wavelength, altitude, and meteorological range required for Figure 4, and 
the twenty tabulations (surface to 5 km) in Tables 3.1 to 3.20. 


7. SUMMARY AND CONCLUDING REMARKS 


To formulate an atmospheric attenuation model with meteorological ranges for 
a haze regime, it was necessary first to define the limits of the haze regime. Fol- 
lowing this, eight meteorological ranges were selected within the regime and the 
surface aerosol attenuation coefficient distribution in wavelength was determined 

| for each meteorological range. From the surface aerosol attenuation coefficients, 

the By (nA, V,) then were computed for km intervals to 5 km altitude (typical depth 

of the aerosol mixing layer) by applying a scale height derived for each meteoro- 

| logical range, Optical thickness values required for the model were computed with 

| Eq. (11), Finally, parameters from an earlier published attenuation model (Elter~ 
man, 1968) were combined with those derived in this paper, in order to provide 
continuity to 50 km altitude. 

The shortest wavelength used in this model ig 0.27 microns. The use of shorter 
wavelengths would required the treatment of O, absorption and its attendant uncer- 
tainties, ‘ihe longest wavelength used is 2.17 microns, Calculations for longer 
wavelengths are complicated by the presence of absorption bands of H,9, CO, and 
their wings, Also, at longer wavelengths and low altitude haze conditions, absorp- 
tion by the aeroso) itself is an unknown factor. In between, a total of 20 wavelengths 
is chosen, within the atmospheric windows and for the ultraviolet region where 
ozone absorption is important (Table 2), If required, a satisfactory interpolation 
for the optical parameters can be made between wavelengths in the region 0,27 to 
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Combining Eqs. (8) and (9) and integrating, 
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about ly because light extinction in this spectral region is caused primarily by 
scattering and ozone absorption and both processes are slowly varying functions of 
the wavelength, This is true of the extinction coefficients, B ext (*)» as well as their 
components, By (A), 6. (A), and generally Bg (A). 

Beyond 1p, the computations did not include molecular absorption. Therefore, 
the tables for wavelengths 1.06, 1.26, 1.67, and 2.17. represent the IR windows 
only, The presence of absorption bands due to H,O0 and other gases does not permit 
interpolation between 1.06 < As 2.17, the near IR region considered in this model, 
unless the interpolation is limited to the Rayleigh and aerosol parameters (no ozone 
absorption present). 


8. TABULATIONS 


The tabulations that follow are in computer notation, For example, read 
5.96 - 2 = 5.96X 1072 and 5,96 2 = 5.96 X 107, 

The format deals systematically with a multiplicity of variables, thus permit- 
ting a variety of exploratory calculations, As an example, the extinction coefficients 
can be used for exploratory transmission calculations, The atmospheric extinction 


coefficient is the sum of all the attenuating components: 


Boxp (sr V,) = BA(h,A) + Bg(h,A) +B (h,A,V,) . (12) 


For horizontal transmission over a path length (d) at selected altitude, wave- 
length, and meteorological range 


Ty(h.A,V,) = exp EP ove (h,A,V,) sd). (13) 


For vertical and slant path transmission from sea level to an altitude of in- 
terest at zenith angle 4, 


To on(h.r ¥) = exp [= Te xcp (FA V,) * pec). (14) 


For vertical and slant path transmission between two altitudes (hy and ho) 
above sea level, 


Tgy(h.A,¥,) " exp “(Te xq(g AV) = Text (yA. V,)) sec 6. (15) 
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For a vertical and slant path transmission from a selected altitude out into 


space, 


Ty-olsAVV,) = exp l-7y,4(h.r, V,) sec 0] . (16) 


When used individually, Rayleigh, aerosol, and ozone parameters are formu- 


lated similarly, 


| 
| 
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